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ABSTRACT 

The LINER galaxy NGC 2639 contains a water vapor megamaser, suggesting the 
presence of an edge-on nuclear accretion disk or torus. This galaxy is thus a good 
candidate for revealing absorption by the torus of any compact nuclear continuum 
emission. In this paper, we report VLBA radio maps at three frequencies and an 
ASCA X-ray spectrum obtained to search for free-free and photoelectric absorptions, 
respectively. The radio observations reveal a compact (< 0.2 pc) nuclear source with 
a spectrum that turns over sharply near 5 GHz. This turnover may reflect either 
synchrotron self-absorption or free-free absorption. The galaxy is detected by ASCA 
with an observed luminosity of 1.4 x 10^^ erg s~^ in the 0.6 - 10 keV band. The X-ray 
spectrum shows emission in excess of a power-law model at energies greater than 4 
keV; we interpret this excess as compact, nuclear, hard X-ray emission with the lower 
energies photoelectrically absorbed by an equivalent hydrogen column of ~ 5 x 10^^ 

_2 

cm . 

If we assume that the turnover in the radio spectrum is caused by free-free 
absorption and that both the free-free and photoelectric absorptions are produced 
by the same gaseous component, the ratio / nidi/ J nudl may be determined. If 
the masing molecular gas is responsible for both absorptions, the required ionization 
fraction is ^ 1.3 x 10~^, which is comparable to the theoretical upper limit derived by 
Neufeld, Maloney & Conger (1994) for X-ray heated molecular gas. The two values 
may be reconciled if the molecular gas is very dense - ^ 10^ cm^'^. The measured 
ionization fraction is also consistent with the idea that both absorptions occur in a hot 
(~ 6,000K), weakly ionized (ionization fraction a few times 10~^) atomic region that 
may co-exist with the warm molecular gas. If this is the case, the absorbing gas is ~ 1 
pc from the nucleus. We rule out the possibility that both absorptions occur in a fully 
ionized gas near lO^K. If our line of sight passes through more than one phase, the 
atomic gas probably dominates the free-free absorption while the molecular gas may 
dominate the photoelectric absorption. 

Subject headings: accretion disks - galaxies: active - galaxies: individual (NGC 2639) 
- galaxies: nuclei - galaxies: Seyfert - radio continuum: galaxies 
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1. Introduction 

Over the last decade, it has become clear that many, perhaps all, active galactic nuclei, are 
surrounded by dusty accretion disks or tori on the pc or sub-pc scale (e.g. Antonucci 1993). When 
viewed equatorially, these tori hide emission from the central regions and they are believed to 
be responsible for the observational difference between broad-line objects (Seyfert Is, broad-line 
radio galaxies and quasars), in which the torus is viewed near polc-on, and narrow-line objects 
(Seyfert 2s, narrow- line radio galaxies), in which the torus is viewed near edge-on. This model is 
supported by a wide range of observational results on narrow-line objects - polarized broad lines 
(e.g. Tran 1995), reddened nuclei at optical wavelengths (Mulchaey, Wilson k. Tsvetanov 1996), 
broad infrared recombination lines (e.g. Goodrich, Veilleux k. Hill 1994), large gas columns to 
the nuclei as inferred from photoelectric absorption of soft X-rays (e.g. Turner et al. 1997) and 
bi-cones of ionized gas aligned with the radio ejecta (e.g. Wilson k, Tsvetanov 1994). 

The obscuring material, which appears to take the form of gas clouds in a geometrically 
thick torus or a warped, thin disk, is illuminated by the central UV and X-ray continuum source. 
Calculations of the physical and chemical properties of the gas have been made by Krolik & 
Lepp (1989) and Neufeld, Maloney & Conger (1994, hereafter NMC). The UV and soft X-rays 
are expected to be absorbed in thin layers at the surface of the clouds, but hard X-rays may 
penetrate and heat the interiors. This X-ray heated gas may possess a two phase structure, in 
which an atomic phase at T ~ 5,000 - 8,000K coexists with a molecular phase at T ~ 600 2,500 
K (e.g. NMC). The ionization fraction is expected to be ^ 10^^ in the molecular region and a 
few times 10~^ in the atomic region, but the exact values are sensitive to the X-ray flux and gas 
density and to details of the chemical and physical processes in the gas. Clearly, an observational 
determination of the ionization fraction would be of value and is the goal of the present paper. 

The ionization fraction may be estimated, in principle, by comparing the emission measure, 
/ nidi (from measurements of free-free absorption of the nuclear radio emission) with the total 
equivalent hydrogen column, / nndl (from measurements of photoelectric absorption of nuclear 
X-rays), both through the disk or torus. To provide maximum path length for these absorptions, 
the disk needs to be oriented close to edge-on. Studies of H2O megamascr emission from the 
nucleus of NGC 4258 reveal that the maser emission arises in a thin Keplerian accretion disk which 
is viewed very close to edge-on (Watson & Wallin 1994; Miyoshi et al. 1995). More recent VLBI 
mapping of other H2O megamasers shows that, in almost all cases, the maser emission traces 
a line on the sky with kinematics consistent with an edge-on, rotating disk (Greenhill & Gwinn 
1997; Greenhih, Moran & Herrnstein 1997; Trotter et al. 1998). Thus galaxies with detected H2O 
megamaser emission provide the best opportunity for measuring the ionization fraction of the 
circumnuclear disk or torus. 

NGC 2639 has a LINER-type nucleus (e.g. Ho, Filippenko & Sargent 1993) and H2O 
megamaser emission (Braatz, Wilson k Henkel 1994). Although a VLBA map of the maser 
emission is not yet available, the systemic maser emission has been found to drift redward at 



- 4 - 



a similar rate to that seen in NGC 4258 (Wilson, Braatz & Henkel 1995). In NGC 4258, this 
redward drift is known to result from the centripetal acceleration of clumps of masing gas on 
the near side of the edge-on accretion disk, as the gas passes in front of the nuclear radio source 
(Herrnstein et al. 1997). By analogy with NGC 4258, it may be argued that the redward drift in 
NGC 2639 arises in the same way and that it too contains an accretion disk viewed very close to 
edge-on. We therefore chose NGC 2639 for a search for absorption by the putative disk. 

Throughout this paper, we adopt a velocity of NGC 2639 with respect to the microwave 
background radiation of Ysk = 3,434 km (de Vaucouleurs et al. 1991) and a Hubble constant 
of 75 km s^^ Mpc^^, giving a distance of 45.8 Mpc and a scale of 222 pc per arc second. 

2. Observations 

2.1. VLBA Radio Observations 

NGC 2639 was observed as part of our survey of Seyfert galaxies on May 31, 1996, using the 
ten antennas of the VLBA^ at 1.7, 5.0, and 15 GHz. Baseline lengths were from 130 to 5000 km, 
and NGC 2639 was observed for 57 min at 1.7 GHz, 48 min at 5.0 GHz, and 46 min at 15 GHz. 
We recorded left circular polarization with 32-MHz bandwidth and 2-bit sampling. Calibration 
and imaging were done with the AIPS software, using standard methods. The flux density scale 
was calibrated using standard VLBA antenna gains and measurements of T^ys made every 1-2 
minutes. 

We phase referenced to the nearby source J0832-I-4913 (2.1° away, adopted J2000 position RA 
= 08^^ 32"^ 23^21671, Dec = 49° 13' 21". 0388 ±0.5 mas [Eubanks 1995, private communication]), 
with cycle times of 9 min at 1.7 GHz and 5.0 GHz, and 2 min at 15 GHz. After phase referencing, 
the data were imaged, deconvolved, and self calibrated, with convergence achieved after two 
iterations of phase-only self calibration. This yielded high-resolution images (from untapered u-v 
data) with FWHM beam sizes of 6.4 mas x 4.6 mas in p. a. 159° at 1.7 GHz, 2.0 mas x 1.7 mas 
in p. a. 169° at 5 GHz, and 1.2 mas x 0.75 mas in p. a. 105° at 15 GHz. The r.m.s. noise on 
these images is n = 0.14, 0.35 and 0.26 mJy (beam area)^^ at 1.4, 5.0 and 15 GHz, respectively, 
which is similar to the expected thermal noise at 1.4 and 15 GHz, but somewhat higher than the 
expected value (0.15 mJy (beam area)~^) at 5.0 GHz. For measuring accurate spectral indices, 
we tapered the array at 5 and 15 GHz to produce approximately the same beamwidth as for the 
untapered 1.7-GHz array, using spacings from 0.5, 1.2 and 4.2 MA to 40 MA at 1.7, 5, and 15 
GHz, respectively. A restoring beam of 6.4 mas x 4.6 mas (FWHM) was then used at all three 
frequencies. The rms noise levels in the tapered images were 0.69 mJy (beam area)~^ at 5 GHz 
and 0.77 mJy (beam area)~^ at 15 GHz. 



^The VLBA is part of the National Radio Astronomy Observatory, a facility of the National Science Foundation, 
operated under cooperative agreement by Associated Universities, Inc 
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The flux-density measurements have uncertainties due to flux-scale calibration, thermal noise, 
and deconvolution errors. The flux scale has a nominal uncertainty of 5% at each frequency 
(uncertainties quoted throughout this paper arc r.m.s.). Deconvolution effects were estimated 
approximately by experimenting with different depths of cleaning and sizes of the clean boxes. 
The test was performed on NGC 1068, which was observed during the same run and in the same 
way as NGC 2639. A grid of nine different sets of reasonable parameters produced a distribution 
of total flux densities that had a standard deviation of 6.8%, and we adopt this estimate as the 
error induced by deconvolution. Resolution effects should be minimal in the measurement of 
spectral indices because we have matched the beamwidths at all frequencies and because the 
source is only slightly resolved. Combining all the above effects in quadrature, the uncertainty on 
the flux-density measurements is ^{O.OSS)'^ + n'^ mJy, where S is the flux density and n is the 
noise given above for each frequency. 

2.2. ASCA X-ray Observations 

NGC 2639 was observed by ASCA (Tanaka, Inoue & Holt 1994) on 16 April, 1997. ASCA 

carries two pairs of instruments known as the solid-state imaging spectrometers (SIS; hereafter 
SO and SI) and the gas-imaging spectrometers (CIS; hereafter G2 and G3). The SIS data 
were obtained in 1-CCD mode, converted to BRIGHT mode, and hot and flickering pixels were 
removed. Good time intervals for the SIS and GIS were determined following Weaver et al. (1994), 
with additional SIS screening to remove times <60 s after South Atlantic Anomaly passages and 
transitions from satellite day to satellite night. 

Data were extracted from within circular regions centered on the source with radii 2.3' and 
5.3' in the SIS and GIS, respectively (larger source regions for the SIS do not improve the signal). 

The SIS background was obtained from an annular ring of width ~ 1.5' on the source chip while 
the GIS background was obtained from circular regions in the source field chosen to avoid a nearby 
serendipitous X-ray source. NGC 2639 has a total ASCA count rate of 0.01 to 0.02 counts s^^ 
depending on the detector, 70% to 81% of which is made up of background counts. Subtracting 
this background yields maximum count rates of 0.01 and 0.006 counts s~^ (keV)~^ for SO and 
G3 respectively. We measure smaller count rates in SI and G2, in which NGC 2639 is positioned 
about 8' from the optical axis. The effects of vignetting cause the galaxy to be indistinguishable 
from the background in these detectors above 3 keV. In SO and G3, NGC 2639 is ~ 5' from the 
optical axis and thus less affected by vignetting. Therefore we consider only data from SO and G3 
in our spectral fits. 

For a weak source like NGC 2639, incorrect background subtraction can cause spurious results. 
Therefore we have examined background taken from blank-sky fields in addition to background 
taken from near the source (described above) . The choice of background makes little difference for 
SO, but has a significant effect on the data from G3. We therefore performed two sets of spectral 
fits with the different G3 backgrounds. The qualitative features of the spectrum do not depend on 
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which background is used and the errors quoted (Section 3.2 and Table 2) include the uncertainty 
in the background. 

3. Results 

3.1. Radio Structure and Spectrum 

Radio fluxes of the nuclear source and beam sizes are given in Table 1 while the radio 
spectrum is shown in Fig. 1. Our 5 GHz flux of 47 it 4 mJy is higher than both a previous VLBI 
measurement of 27 it 4 mJy (Hummel et al. 1982, made between June 1980 and April 1981) and 
the peak brightness of 23.4 it 0.9 mJy beam~^ in a VLA map (Ulvestad & Wilson 1989, made on 
Feb 3 - 4 1985), indicating the source is variable. The source is unresolved at 1.7 and 5.0 GHz, 
but extended at 15 GHz. The peak flux density in our highest resolution 15 GHz map is 33.8 mJy 
(beam area)~^, significantly lower than the total flux of 45 it 4 mJy. The FWHM deconvolved 
source size is 0.70 x < 0.15 mas (0.16 x < 0.03 pc) with major axis in p. a. 111°. For comparison, 
VLA observations reveal an incompletely resolved triple source with overall length 1.4" (310 pc) 
in p. a. 105° (Ulvestad &: Wilson 1989). The agreement between the p.a.'s of the sub-pc scale and 
the hundreds of pc scale radio emission strongly suggests that the extension found with the VLBA 
observations represents the inner part of the jet presumed to fuel the VLA-observed structure. 
The spectrum of the nuclear source (Fig. 1) is flat between 5 and 15 GHz (a^ = +0.04 it 0.08, 
S (X v~'^) and turns over sharply below 5 GHz (af ^ = —1.8 =b 0.1). We now discuss various 
interpretations of the emission mechanism of the nuclear radio source and the cause of the low 
frequency turnover. 

The brightness temperatures, T^, of the nuclear source are > 1.5 x lO^K, > 9.8 x lO^K and 
> 2.9 X lO'^K at 1.7, 5.0 and 15 GHz, respectively. For comparison, the brightness temperature of 
the source SI in NGC 1068, which is believed to coincide with the nucleus and may be thermal 
m origin, is of order a few xlO^ K to a few xlO^ K at 8.4 GHz (Gallimore et al. 1997; Roy et 
al. 1998). The high brightness temperature of the source in NGC 2639 makes a thermal origin 
unlikely. 

The three processes that could, in principle, be responsible for the low frequency turnover are 
synchrotron self-absorption, Razin-Tsytovich suppression and free-free absorption, and we discuss 
each in turn. We can ascribe an effective temperature to electrons of a given energy through 
the equation 3A;Te = ^(.m(,(? or Tg = 2.0 x 10^7eK, where 7e is the relativistic gamma factor of 
the radiating electrons, me is the electron mass, c is the speed of light and k is Boltzmann's 
constant (e.g. Longair 1981). For a synchrotron self-absorbed source, = Tg. The brightness 
temperature measurements are only lower limits and are obviously consistent with a synchrotron 
self-absorption interpretation. The source is, however, resolved in one dimension at 15 GHz. If 
we adopt a 15 GHz size of 0.16 x 0.03 pc (where the minor axis extent is taken to be equal to 
its upper limit) and assume the size is the same at 5 GHz (where the spectrum begins to turn 
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over), then Ti,{5 GHz) ~ 2.9 x 10^°K, which allows the possibility of synchrotron self-absorption. 
The mean magnetic field strength would be ~ 0.5 gauss (using equation (1) of Kcllermann &; 
Pauliny-Toth 1981), which is similar to the equipartition value of ~ 0.3 gauss (calculated assuming 
that the total energy in cosmic rays is twice that in electrons, the radio source is optically thin at 
15 GHz, the electron energy spectrum extends from radiated frequencies 10 MHz to 100 GHz with 
index 7 = 2.5 {N{E) oc E~'^), and that the source is ellipsoidal in shape with axes 0.16 x 0.03 x 
0.03 pc). We should, however, bear in mind that the resolved flux at 15 GHz is only ~ 25% of 
the total flux, so most of the observed emission may come from unresolved regions with higher 
brightness and weaker magnetic fields. 

Razin-Tsytovich suppression occurs at frequencies below ~ 20ne/B±, where Ue is the thermal 
electron density (in cm^'^) within the synchrotron-emitting region and i?_L is the perpendicular 
component of the magnetic field (in gauss) . Comparison with the equation for free- free absorption 
(see below) shows that the Razin-Tsytovich cut-off frequency will lie above the free-free absorption 
frequency only if B± < 4 x 10~^L~^/^ gauss, where L is the dimension of the absorbing region in 
pc (e.g. Moffet 1975). We can take the measured maximum extent of the source at 15 GHz - 0.16 
pc - as an upper limit for L, so for Razin-Tsytovich suppression to occur at a higher frequency 
than free-free absorption, the magnetic field would need to be B± < 1 x 10^^ gauss, which is 
more than three orders of magnitude below equipartition. We conclude that Razin-Tsytovich 
suppression is an implausible explanation for the low-frequency turnover in NGC 2639. 

To investigate the possibility that free-free absorption by ionized gas along the line of sight 
to the radio source is responsible for the low frequency cut-off, we assume the observed spectrum 
may be described by: 

S{u) = aiz-'^e-^M (1) 

with 

t{u) = 8.235 X IQ-'^T-^-^^u-'^-^ J nidi (2) 

where a is a constant, r(i/) is the optical depth to free-free absorption, (K) and Ug (cm^^) are 
the electron temperature and density of the absorbing gas, u is the frequency (GHz) and I is the 
path length in pc (e.g. Mezger Sz Henderson 1967). The three observed flux densities may be 
used to solve for the three unknowns (a, a and r) in equation (1), given the known dependence of 
r on in equation (2). The result for the spectral index and optical depth is a = 0.25 and r(5 
GHz) = 0.25. The value of a should not be taken too seriously (a flux measurement at a higher 
frequency is needed to tie down the optically thin part of the spectrum), but r is well defined by 
the observed spectrum. The spectrum from this model is shown by the line in Fig. 1. The derived 
value of r(i/) implies 

J nidi = 2.3 X lO^rJ/^cm-^pc (3) 
where Te4 = Te/W^ K. We discuss the nature of the absorbing gas in Section 4. 
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3.2. X-ray Spectrum 

We first modelled the X-ray spectrum as a single, uniformly-absorbed power-law, but found 
the fit to be poor (model A, Table 2), with the data showing an excess over the model above 4 
keV. This excess emission is independent of the choice of G3 background and probably represents 
a heavily absorbed, hard component. A model (model B) comprising a single power-law spectrum 
absorbed by an intrinsic medium covering 93% of the source (i.e. "partial covering") and the 
Galactic column density (Nh = 3 x lO^'^ cm~^) covering the whole source, provides an excellent 
description of the data (again, independent of the choice of G3 background). The part of the X-ray 
source covered by only the Galactic column could be scattered or intrinsically extended X-ray 
emission. The best estimate of the photon index is F = 2.4 it 0.4. Constraining the power-law 
index of the photon spectrum to F = 1.9 (close to the 'canonical' value for Scyfcrt galaxies) does 
not significantly worsen the fit (model C). We also tried a model that consists of a power law 
absorbed by the Galactic column plus a Gaussian function to represent the Fe Ka emission often 
seen in active galaxies (model D). Such a model describes the high-energy excess as well as models 
B and C; however, the line energy is much too low to be Fe Ka emission at the redshift of NGC 
2639 (Table 2). 

Both models B and C require an intrinsic column density / nndl = (4 — 5) x 10^^ cm~^, 
although columns half to almost three times as large are acceptable at the 90% confidence level 
(Table 2) . The observed SO and G3 spectra and model B folded through the instrumental response 
are shown in Figure 2. 



4. The Ionization Fraction 

A number of assumptions must be made in calculating the ionization fraction. First, we 
assume that the low frequency turnover in the radio spectrum results from free-free absorption. 
There is now strong evidence that the VLBI-scale radio sources in several Seyfert galaxies suffer 
free-free absorption at GHz frequencies (Halkides et al. 1998; Roy et al. 1998; Ulvestad et al. 
1998), and our spectrum of NGC 2639 is consistent with this process (synchrotron self-absorption 
is, however, a viable alternative - see Section 3.1). Second, the lines of sight to the hard X-ray 
source and the compact radio source must be the same. While the resolution of the VLBA is high 
and the radio source is compact (< 0.2 pc), coincidence of the two sight lines is not guaranteed. 
For example, if some of the radio emission originates from a compact jet, it may be spatially 
separated from the presumed nuclear hard X-ray source. Third, we assume that the same gas 
is responsible for both free-free and photoelectric absorptions. This is not necessarily so, even if 
the two lines of sight coincide: it is possible that the photoelectric absorption is dominated by 
weakly ionized atomic or molecular gas which is dense or far from the central ionizing source, 
while the free-free absorption may result from highly ionized gas of low density or close to this 
source. The apparent ionization fraction would then represent a suitably weighted average along 
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the line of sight and is not necessarily representative of this parameter in a particular gaseous 
component. Later in this Section (Section 4.4), we explore the consequences of abandoning this 
third assumption. 

Our measurements imply / n^dl/ J nfjdl = 140 T^4^^ cm~^, with a formal uncertainty of a 
factor of 2 (dominated by the uncertainty in / nndl). For a uniform absorbing medium, the 
ionization fraction is then: 

Xe = Ue/UH = 1.2 X lO-^rO/75^-0.5 (4) 

where nnio = nn/W^^ cm~^. The uncertainty in / nndl translates into an uncertainty of ~ 1.4 in 
Xe- Within the context of a model in which dense gas is heated and photoionized by a nuclear UV 
and X-ray source, the absorption could occur in a warm, weakly ionized molecular phase, a hot, 
weakly ionized atomic phase or a hot, highly ionized atomic phase (listed in order of decreasing 
distance from the ionizing source). We consider these three possibilities in turn. 

4.1. Warm Molecular Phase 

For water vapor maser emission, we require (e.g. Elitzur, Hollenbach & McKee 1989) 
^ lO^'^ cm~^ (to avoid collisional quenching) and gas temperatures Tg ~ 600K (for coUisional 
pumping). Thus, if the observed column of absorbing gas comprises the masing region, equation 
(4) implies Xg ^ 1.3 x 10~^ (with the equality for ra^/j = ^^^^ cm~^), assuming Tg = Tg. This 
lower limit is comparable with the upper limit of Xg ^ 10~^ derived theoretically for X-ray heated 
molecular gas by NMC. The theoretical, column-averaged ionization fraction may be as high as 
several x 10~^, depending on the X-ray flux, gas column and pressure (P. R. Maloney, private 
communication). Thus warm molecular gas could potentially provide both the free- free and 
photoelectric absorptions, but only if the density is high (nj^j ^ 10^ cm~^). 

4.2. Hot, Weakly Ionized, Atomic Phase 

If, on the other hand, the absorption occurs in the hot (~ 6,000K) atomic region in pressure 
balance with the molecular region (e.g. NMC), equation (4) implies Xg ^ 2.7 x 10^^, which is 
consistent with the few times 10~^ expected by NMC. If our line of sight passes through this 
atomic region, its free-free absorption should dominate the molecular phase (NMC). Demanding 
NMC's ionization fraction of ~ 2 x 10~^, our data imply uh ~ 2 X 10^ cm~^ in the atomic region 
and a path length of Nh/tih ~ 1 pc for the absorbing gas. 

The intrinsic 1 - 100 keV X-ray luminosity of NGC 2639 is Lx = 1.5 x lO'^^ ^j-g g-i (fQj. 
r = 1.9) or 1.1 x 10^^ erg s^^ (for T = 2.4). The X-ray flux incident on the gas is Fx = lO^Fs erg 
cm~^ s~^, where F5 = 0.084(Lx/10'*^ erg s~^) (i?/pc)~^ and R is the distance from the X-ray 
source to the illuminated surface of the gas. For the masing molecular gas and the absorbing 
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atomic phase to co-exist, NMC's Fig. 1 implies F^/{N2ipii) — 10, where Nh = lO^^A/24 cm~^ 
and p/k = lO^^pn K cm^^ {p is the pressure). Taking N2A = 0.5 (from the X-ray absorption), 
uh = 2 X 10^ cm~^ (see above) and T = 6, COOK for the atomic region, wc find F5 ~ 0.06 and 
R ~ 1.2 — 1.4 pc. This size scale corresponds to ~ 5 — 6 mas, which will be easily resolvable in 
planned VLB A observations of the H2O megamaser. 

4.3. Hot, Fully Ionized, Atomic Phase 

Next we investigate the hypothesis that a fully ionized {x^ — 1) gas at lO^K is responsible for 
both absorptions. In this case, equation (4) gives Ue ~ riH ~ 140 cm^^ and the absorber path 
length I = Nh/uh ~ 1 kpc. For a Stromgren sphere geometry, the required number of ionizing 
photons is 

iV* = aBnlAirf/?, (5) 

{aB is the total hydrogen recombination coefficient in case B), giving A^* = 9.6 X 10^6 photons 
s~^. This rate of ionizing photons may be compared with = 3.3 x 10^^ photons s~^ obtained 
by extrapolating the measured hard X-ray flux down to the Lyman limit with F = 2.4, the best 
estimate of the photon index of the hard X-ray source. This value of F is also typical of radio-quiet 
quasars between the Lyman limit and 1 keV (e.g. Laor et al. 1997); however, the spectra of 
LINERs in the ionizing UV arc poorly known, so our estimate is very uncertain. Nevertheless, 
given that the required value of N^, is more than three orders of magnitude larger than this 
estimated value, we rule out the possibility that both free-free and photoelectric absorptions occur 
in the fully ionized phase. 

4.4. Multiple Phases 

Lastly, we abandon the assumption that the same gas is responsible for both types of 
absorption and consider the possibility that the free-free absorption occurs in the fully ionized gas 
at lO^K while the photoelectric absorption is dominated by another, weakly ionized component. 
Taking N^, = 3.3 x 10^^ photons s~^ from the extrapolation of the X-ray spectrum, T^a = 1 and 
combining equations (3) and (5), we have rig ~ 1,000 cm^^ and I 20 pc, suggestive of a large, 
dense cloud in the narrow line region. The implied column density through the fully ionized layer 
is Nh ~ 7 X 10^^ cm^^, a factor of 7 lower than inferred from the X-ray spectrum, so a weakly 
ionized column of ~ 4x 10^^ cm~^ would also have to be present. The fully ionized gas would 
produce an H/3 flux of 1.2x10"^^ erg cm~^ s~^, which is factor of 8 larger than the observed, 
reddening-corrected H/5 flux from the nucleus of NGC 2639 (Ho, Filippenko & Sargent 1993), so 
most of this gas would have to be obscured at optical wavelengths in this interpretation. 
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4.5. Synopsis of Absorption Results 

In summary, our data are consistent with the absorptions occurring in any of: (i) a warm 
molecular region of density ^ 10^ cm~^ and ionization fraction several x 10~^, (ii) a hot atomic 
region of density ~ 10^ cm~^ and ionization fraction a few x 10^^ which co-exists with the 
masing molecular gas about 1 pc from the miclear X-ray source, or (iii) a combination of atomic 
and molecular regions, with the hot atomic region dominating the free-free absorption and the 
molecular region dominating the photoelectric absorption. 

5. Conclusions 

We have found that the spectrum of the nuclear radio source in NGC 2639 turns over below 
^ 5 GHz and shown that this effect may result from either synchrotron self-absorption or free-free 
absorption. If the latter process is responsible, the implied value of / ri^dl may be combined 
with the equivalent hydrogen column, / nndl (derived from the X-ray spectrum), to provide the 
ionization fraction of the gas in terms of the electron temperature and hydrogen density, assuming 
the same gas is responsible for both kinds of absorption. The lower limit to the ionization fraction 
obtained by assuming the absorptions occur in the masing molecular gas is comparable to the 
upper limit to the ionization fraction derived theoretically by Neufeld, Maloney &: Conger (1994) 
for warm molecular gas heated by a nuclear hard X-ray source, implying that the molecular gas 
must be dense {'^ 10^ cm~^) if it is to provide the free-free absorption. The required ionization 
fraction is consistent with that expected in the atomic phase at 5,000 - 8,000 K which may 
co-exist with the molecular component. If our line of sight passes through both phases, the hot 
atomic phase would likely dominate the free-free absorption while molecular gas could dominate 
the photoelectric absorption. We rule out the hypothesis that both types of absorption occur in 
a fully ionized gas at 10^ K. Planned VLBA imaging of the H2O megamaser will help distinguish 
between these various possibilities. 
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Table 1. Flux Densities, Beam Sizes (FWHM) and Position from VLBA Observations 



Frequency (GHz) 





1.7 


5.0 


15.3 


Flux Density (mJy) 


6.9±0.6 


47±4 


40±3 


Beam Size (mas) 


6.4x4.6 


2.0x1.7 


1.2x0.75 


Flux Density (mJy) 




47±4 


45±4 


Beam Size (mas) 




6.4x4.6 


6.4x4.6 



J2000 position of peak radio flux: RA = OS'' 43"^ 
38^07788, Dec = 50° 12' 20".0044 ±0.9 mas (absolute 
position uncertainty). 
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Table 2. Models of the X-ray spectrum 



Model" 






N//(soft) 

( 10^2 cm-2) 


N//(liard) 
( 10^2 cm-2) 


xVd.o.f. 
SB^ 


xVd.o.f. 
BB^ 


A 


1.8 (1.4 - 


2Af 


(0 - 0.19)^* 




56.5/45 


55.3/45 


B" 


2.4 (1.9 - 


2.8f 


omf'9 


42 (19 - 98)'^ 


44.2/44 


45.6/44 


C 


1.99 




0.03-^'9 


53 (22 - 130)<^ 


47.7/45 


48.5/45 




1.99 




0.03-^'9 




46.7/45 


42.0/45 



'^Models arc A: uniformly-absorbed power law; B: power-law source partially covered by N//(hard) 
and fully covered by N//(soft); C: As model B, but with T fixed at 1.9; D: power law plus Gaussian 
(intended to represent Fe Ka emission) with Galactic absorption only. 

^Photon index. 

'^SB and BB represent fits with the G3 background taken from the source and blank-sky fields, 
respectively. 

'^The values given represent the 90% confidence ranges for the two free parameters T and Nh 
in models A and B, and the 90% confidence range for the single free parameter Nh in model C. 
Errors on F and N/f include the uncertainty in the G3 background. 

^Cloud covering factor is 0.93 (90% confidence range is 0.83 - 0.98). 
^Galactic hydrogen column density toward NGC 2639. 

^Fixed parameter. 

'^The energy and equivalent width of the Gaussian are 5.5 (90% confidence range is 5.0 - 5.7) 
keV and ~ 3 keV, respectively. The line width is fixed at a = 0.01 keV. 

Note: the 0.6 — 10 keV and 2 — 10 keV observed (absorbed) fiuxes are ~ 6 x 10~^^ and ~ 4 x 
10~^^ erg cm~2 s~^, respectively. The 0.6 — 10 keV and 2 — 10 keV intrinsic (unabsorbed) fiuxes 
are 2± 4 x lO""^^ and ~ 2 x lO""^^ erg cm~2 s~^, respectively. 
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Fig. 1. — The radio spectrum of the nuclear source in NGC 2639. The points are the VLB A 
measurements, with the radii of the circles equal to the r.m.s. uncertainties in flux density. The 
line represents a model comprising a source with a power-law spectrum which is free-free absorbed 
by intervening ionized gas (see Section 3.1). 

Fig. 2. — The observed ASCA spectra of NGC 2639 from detectors SO and G3 along with model B 
(see Table 2) folded through the instrumental response. Thick lines are SO data and model; thin 
lines are G3 data and model. 



